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ABSTRACT 


An  analytical  study,  using  the  STAGSC-1  computer  code,  was  conducted 
on  a  graphite/epoxy  (AS4/3501  )  composite  cylindrical  panel  acting  under- 
compressive  loading  considering  cutouts  positioned  at.  critical  locations 
within  the  panels  surface  area.  The  study  also  investigated  the  material 
degradation  effects  of  temperature  and  moisture  on  the  collapse  load  of 
these  cylindrical  panels.  Two  temperatures  (80  and  250  °F),  in  both  a 
symmetric  and  unsymmetric  moisture  condition,  were  investigated.  The 
overall  material  degradation  characteristics  were  investigated  by 
degrading  the  E2  and  G12  moduli  of  the  individual  plies. 

It  was  found  that  the  effects  of  moisture  and  temperature  (material 
degradation  effects)  can  greatly  reduce  the  collapse  load  of  the  cutout 
panel.  A  15X  decrease  in  the  collapse  load  occurred  at  the  saturated 
moisture  conditions  and  elevated  temperatures.  The  material  degradation 
effects  were  found  to  produce  the  same  results  no  matter  where  the  cutout 
was  located.  The  panel  with  cutouts  was  less  effected,  at  the  saturated 
moisture  condition  and  elevated  temperature,  then  that  of  a  panel  with  no 
cutout. 

It  was  also  found  that  the  collapse  load,  as  expected,  was  dependent 
upon  the  location  of  the  cutout  within  the  panel.  There  was  a  decrease  in 
the  collapse  load  by  as  much  as  26.2%  when  the  cutout  was  moved  from  the 
center  of  the  panel  to  a  position  closer  to  the  side  edge.  There  was 
very  little  change,  IX,  when  the  cutout  was  positioned  closer  to  the 
loading  edge  of  the  panel.  Ply  layup  was  also  of  great  importance.  There 
was  a  21. 3X  increase  or  decrease  iri  the  collapse  load  if  a  certain  ply 
layup  (0/±45/90)s,  (symmetrical  ly  exchanging  of  the  second  and  third 
plies)  was  chosen  for  the  cutout  positioned  near  the  side  supports. 


I.  INTRODUCTION 


Background 

Large,  lightweight  shell  structures  are  used  extensively  in  aircraft 
for  both  strength  and  stability.  These  shell  structures  are  greatly 
influenced  by  either  reinforced  or  nonreinforced  cutouts  and  material 
degradation  conditions.  Due  to  these  influences  and  with  the  recent 
advances,  a  vast  number  of  these  shells  are  made  from  composite 
materials.  The  composite  materials  are  desired  because  of  the  high 
strength  to  weight  ratio  when  compared  to  conventional  materials. 
Therefore,  from  a  practical  point  of  view,  there  is  a  need  for  further 
research  into  composite  shell  structures  and  the  effect  of  cutouts  and 
material  degradation  conditions  on  the  stability  of  the  shell  structures. 

Very  few  references  can  be  found  dealing  with  buckling  of  composite 
panels  and  plates  under  axial  compression  with  cutouts.  In  contrast, 
numerous  studies  have  been  carried  out  on  isotropic  plates  and 
cylindrical  panels  without  cutouts  [1-6].  These  studies  centered  around 
various  boundary  conditions  and  their  affect  on  the  buckling  load.  If 
one  considered  the  area  of  axially  loaded  isotropic  panels  with  cutouts 
[7-11],  mainly,  centered  circular  cutouts  were  analyzed  numerically  by 
finite  difference  techniques  to  determine  the  buckling  and  bifurcation 
loads . 

As  mentioned  earlier,  composite  panels  have  become  an  important  part 
of  aerospace  structures  due  to  their  excellent  strength  to  weight  ratio. 
Composite  cylindrical  panels  without  cutouts  have  been  studied  [12-16], 


but  very  little  work  can  be  found  dealing  with  composite  panels  with 
cutouts.  Janisse  [17]  and  Lee  [13]  are  two  authors  who  have  studied  the 
problem  of  cutouts  in  composite  cylindrical  panels.  Janisse  [17]  found 
that  the  collapse  characteristics  of  composite  cylindrical  panels  are 
dependent  on  ply  layup  and  size  of  the  cutout.  Lee  [18],  in  a  follow-up 
investigation,  studied  the  effects  of  the  cutout  aspect  ratio  and 
concluded  that  as  the  surface  area  of  the  cutout  increased  the  buckling 
load  decreased  with  the  effects  of  nonlinearity  becoming  more  pronounced. 
These  are  the  only  studies  which  deal  specifically  with  composite 
cylindrical  panels  containing  cutouts.  This  shows  that  very  little  is 
known  about  the  effect  of  cutouts  and  further  study  into  this  area  is 
needed . 

Temperature  and  moisture,  the  main  emphasis  of  an  environmental 
study,  is  another  important  area  of  interest  in  composite  shell 
structures.  Snead  [19],  while  examining  the  effects  of  moisture  and 
temperature  on  the  instability  of  cylindrical  composite  panels,  concluded 
that  there  was  a  definite  degradation  of  the  bifurcation  load  due  to 
elevated  temperatures  and  moisture  exposure.  Various  other  authors 
[20-23]  have  looked  at  the  environmental  effects  and  have  come  to  the 
same  conclusion.  There  is,  however,  a  lack  of  research  into  how  elevated 
temperature  and  moisture  exposure  (material  degradation  effects)  relates 
to  the  collapse  of  a  composite  cylindrical  panel  with  a  cutout. 

With  a  need  for  faster,  lighter  aircraft  and  the  increasing  interest 
in  space,  a  better  understanding  of  composite  cylindrical  shells  with 
cutouts  is  needed.  All  of  the  previous  cutout  studies  that  were  found 
dealt  with  centrally  located  cutouts  and  standard  conditions.  This 


thesis  addresses  the  issue  of  eccentrically  located  cutouts  in  composite 


cylindrical  panels,  and  the  effect  of  material  degradation  due  to 


temperature  and  moisture  on  the  collapse  load  of  an  axially  loaded  panel. 


Ob jec  ti ve 


The  ira-jor  purpose  of  this  thesis  is  to  study  the  effects  of  material 


degradation  on  panels  with  cutouts.  To  observe  this  effect,  one  first 


has  to  establish  the  critical  positions  of  specific  size  cutouts  within 


the  panel.  In  order  to  do  this,  the  effect  of  a  cutout  located  at 


different  positions  within  the  composite  cylindrical  panel  subjected  to 


an  axial  loading  will  be  studied  analytically  using  the  Structural 


Analysis  of  General  Shells  (STAGSC-1)  computer  program.  Then  secondly, 


material  degradation  effects  due  to  temperature  and  moisture  at  the 


critical  cutout  locations  will  be  analyzed.  Upon  completion,  this  thesis 


will  allow  for  a  better  understanding  of  axially  loaded  composite 


cylindrical  shells  with  cutouts  under  less  then  perfect  conditions. 
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Approach 


Due  to  the  presence  of  the  geometric  discontinuity  (cutout),  a 


nonlinear  analysis  was  used.  This  nonlinear  analysis  is  part,  of  the 
5TAGSC-L  finite  element  computer  program  [26-27],  developed  by  Lockheed, 
which  incorporates  a  modified  Newton-Raphson  iteration  technique  along 
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with  large  displacement  theory  and  moderate  rotations  to  study  the 
collapse  of  cylindrical  panels  with  cutouts.  The  cylindrical  shell 
geometry  used  by  the  STAGSC-1  computer  program  is  shown  in  Fig  1.1. 

Eccentrically  located  cutouts  were  investigated  by  positioning  the 
cutout  at  four  different  locations  as  shown  in  Fig.  1.2.  Both  a  2"  x  2 
cutout  (2.732  of  panels  total  area)  and  a  4“  x  4”  cutout  (11.11)  were 
examined  along  with  the  effect  of  interchanging  ply  positions. 

Material  degradation  effects  on  a  cutout  located  at  different 
positions  within  the  panel  were  also  investigated.  Cutout  locations  A 
and  C,  Fig.  1.2,  were  chosen  for  this  analysis.  For  the  cylindrical 
panel,  a  symmetric  and  a  unsymmetric  moisture  concentration  condition  w 
used.  Two  time  units  were  chosen  for  each  of  the  moisture  conditions;  a 
steady  state  distribution  and  a  near  zero  moisture  absorption.  These 
time  units  represent  the  two  extremes  of  moisture  conditions.  Finally, 
two  different  temperatures  were  chosen,  80  °F  (300  °K)  and  250  °F  (394 
° K ) .  The  runs  produced  by  the  combination  of  the  fore-mentioned 
conditions  represents  a  near  standard  zero  condition  (room  temperature 
and  near  zero  moisture)  and  that  of  a  high  temperature  and  fully 
saturated  moisture  condition. 


THEORY 


Classical  Lamination  Theory 


An  in-depth  understanding  of  composite  laminae  behavior  is  needed  if 
one  is  going  to  analyze  composite  structures.  Presented  here,  is  a  brief 
overview  of  the  basic  principles.  The  reader  should  refer  to  Ref.  [24] 
and  [25]  for  the  in-depth  development  of  these  relations  and  a  better 
understanding  of  Classical  Lamination  Theory. 

For  a  lamina  of  orthotropic  material,  the  plane  stress  constitutive 
relations  in  principal  material  coordinates  are 
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where  Qij’s  are  the  reduced  stiffnesses  in  terms  of  the  engineering 
constants  and  e  i  and  ^2  are  normal  strains  and  T 1 2  is  the  shear  strain 
for  the  lamina  principle  axes  as  shown  in  Fig.  2.1.  Poisson’s  ratio, 

^ i j ,  and  shear  modulus,  G12,  for  the  1-2  plane  are  defined  as 
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So  far,  orthotropic  materials  have  been  dealt  with  in  their 
principal  material  directions.  Expanding  these  relationships  to  cover 
any  coordinate  system,  with  the  fiber  axis  oriented  at  some  angle  0, 
(Fig.  2.1),  the  stresses  become 
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where 

Mil  =  Miicos4  0  +  2(Qi2  +  2Q66 ) sin2  B  cos2  0  +  U22Sin4  $ 

Q12  =  (Qii  +  M22  -  4Q$s  )sin2  B  cos2  B  +  Qizlsin4  B  +  cos4  B  ) 

M 2 2  -  Qusin4  B  *  2(Qi2  +  2Q66)sin20  cos2  B  +  Q22COS4  6 

Qit  -  (Q11  -  Miz  -  2Qt6)sin  B  cos3  B  +  (Q12  -  Q22  +  2Q6t) 

sin3  B  cos  B  (5) 

Uze  -  (Qii  -  Q 1 2  -  2Q66)sin3  B  cos  B  +  ( 1 2  ~  U22  +  2Qte) 
sin  B  cos3  B 

Ut<>  =  (Mil  +  Mz 2  -  2Qi2  -  2U66)sin20  cos2  B  +  Ueetsin4#  +  cos4  d  ) 


The  above  equations  are  derived  for  individual  lamina  and  in  order 
to  define  the  stress  and  strain  variations  through  the  thickness  (Fig. 
2.2),  each  kth  layer  of  the  multilayered  laminate  can  be  written  as 


10 


I* 


{  a  }k  -  [  a  ]k{  «  h 


By  extending  the  stress-strain  relationships  to  the  multilayered  laminate 
two  assumptions  are  made.  First.,  the  laminate  is  presumed  to  consist  of 


vS 


perfectly  bonded  laminae.  These  laminae  are  infinitesimally  thin  with 


the  displacements  continuous  across  the  laminae  so  that  no  lamina  can 
slip  relative  to  another.  Secondly,  the  Kirchhof f-Love  hypothesis 


applies.  This  hypothesis  sates  that  normals  to  the  mid-surface  remain 
plane  and  normal  to  that  surface  after  bending. 

The  middle  surface  strain-displacement  relations  for  a  cylindrical 


laminate  with  moderately  large  rotations  of  tangents  to  the  panel 
reference  are  given  by  Sanders’  kinematic  relations  without  initial 
imperfections  as 


=  u°,K  -  u ,  x  +  1/2  4>l  +  1/2  02 
:  V °  ,  y  S  V  ,  y  +  W/R  +  1/2  0?  ~  1/2  02 
^  V,x  +  U,y  +  0x  0 y 


where  u,  v,  and  w  denote  the  axial,  circumferential,  and  radial 
components  of  the  mid-surface  displacement.  The  0’s  are  the  rotation 
components  and  in  terms  of  displacement  are 
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where  R  is  the  radius  of  curvature  of  the  panel  at  the  mid-surface. 

The  middle  surface  curvatures  are  then  written  as 

*x  -  0x  ,  X 

ACy  :  0y  i  y  ( 9 ) 

2  (fxy  -  2  AC  y  x  -  0  y  ,  x  +  0  x  ,  y  +  0  /  R 

Using  the  Kirchhof f-Love  hypothesis  with  the  known  middle  surface 
strains  and  curvatures,  the  strain-curvature  relationship  for  a  laminate 
becomes 
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If  the  strain-curvature  relationship,  (Eqn.  10),  which  represents 
the  strain  variation  through  the  thickness,  is  substituted  into  the 
stress-strain  relations,  (Eqn.  6),  the  stress  for  any  kth  layer  can  be 
expressed  in  terms  of  the  middle  surface  strains  and  curvatures  as 
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By  integrating  the  stresses  of  each  laminate  layer,  the  resultant 


forces  and  moments  acting  on  a  laminate  become  for  example 


(12) 


mho  c  '!  ^ r. d  Mx  are  the  force  and  moaent  per  unit,  length  of  the  cross 


section  of  the  laainate  as  shown  in  Fig.  2.3. 


STAGSC-1  Theory 


The  Structural  Analysis  of  General  Shells  (STAGSC-1)  is  a  computer 
program  developed  by  Lockheed  Palo  Alto  Research  Laboratory  [26-27],  The 
program  is  an  energy  based  finite  element  program  that  uses  a 
Newton-Raphson  nonlinear  equation  solver.  STAGSC-1  is  a  general  purpose, 
thin  shell,  structural  analysis  program  used  to  analyze  general  shells 
under  various  static,  thermal,  and  mechanical  loadings.  This  nonlinear 
analysis,  used  to  study  the  collapse  of  shells  with  cutouts,  allows  for 
the  determination  of  the  critical  load. 

According  to  the  energy  method,  a  system’s  total  potential  energy  is 
used  to  derive  its  equilibrium  equations  from  which  stability  can  be 
determined  by  the  solution  to  an  eigenvalue  problem.  A  shell’s  total 
potential  energy  is  equal  to  its  internal  strain  energy  minus  the  product 
of  the  external  forces  and  their  respective  deflections. 

An  element’s  strain  energy  is  given  by  [13]: 


The  expressions  for  the  mid-surface  strains  were  presented  in  Eqn’s.  7 
and  9  and  [N]  was  given  in  Eqn’s.  15  and  16.  In  general,  the  strain 
vector  [  € ] o  is  a  function  of  the  mid-surface  displacements  u,v,w,  the 
first  order  derivatives  of  u,v,w  and  the  second  derivative  of  w. 
Therefore,  we  can  let  a  vector  [d]  represent  this  functional  dependence 
of  [<]o  on  the  displacement  by 

[d]T  =  [u,  U , x ,  U , y ,  V,  V , x ,  V , y ,  W,  W,x,  W,y,  W,xx,  W,xy,  W,yy]  (20) 

By  using  vector  [d],  Bauld  [13]  carried  out  the  integration  of  Eqn.  17 
and  found  that  the  expression  for  strain  energy  is  comprised  of  three 
distinct  parts.  The  first  part,  is  quadratic  in  displacements,  the  second 
cubic  in  displacements,  and  the  third  quartic  in  displacements.  In  a 
finite  element  analysis,  vector  [d]  is  represented  by  an  element’s  shape 
functions  and  nodal  degrees  of  freedom.  Thus,  the  strain  energy  can  be 
expressed  in  terms  of  an  element’s  degrees  of  freedom.  An  element’s 
external  potential  energy  is  obtained  by  using  a  similar  finite  element 
analysis  on  its  external  forces.  The  strain  energy  and  external 
potential  energy  are  combined  to  give  the  total  system’s  potential 
energy,  V,  and  can  be  written  in  the  form  given  by  Bauld  [13]  as: 
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V  3  (1/2  Am  +  1/6  N 1  r  s  +  1/12  N2r5)qrqs  -  Rs qs  (21) 

An  is  the  system’s  linear  stiffness  matrix  with  no  dependence  on  the 
displacement  vector,  q.  Nlr*  and  N2rs  are  matrices  with  linear  and 
quadratic  dependence,  respectively,  on  displacement.  Rs  is  the  surface 
force  vector. 

The  principle  of  total  potential  energy  states  that  the  equilibrium 
configuration  of  a  conservative  mechanical  system  corresponds  to  a 
stationary  value  of  the  total  potential  energy  of  the  system.  Therefore, 
taking  the  first  variation  of  Eqn.  21  and  setting  it  equal  to  zero,  we 
have  a  set  of  nonlinear,  algebraic  equations  of  the  form 


(Ars  +  1/2  N 1 r  j  +  1/3  N2r*)qs  -  R*  =  0 


(22) 


Loss  of  stability  (collapse),  results  when  the  second  variation  of  the 
system’s  total  potential  energy  ceases  to  be  positive  definite,  or 


DET  (Ar  ,  +  Nlr  s  +  N2r  5  )  3  0 


(23) 


Equation  23  is  also  used  by  STAGSC-l  to  solve  the  eigenvalue  problem  of 
the  form 

[A]  +  X [B ]  +  A 2  [ C ]  3  0  (24) 
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where  [B]  and  [C]  represent  nonlinear  stiffness  matrices  in  unknown 


displacements  and  products  of  displacements,  respectively.  For  a  linear 
analysis  the  matrix  C  is  neglected  and  X  is  the  proportionality  constant 
of  a  convenient  load  level  used  in  the  equilibrium  equations  (Eqn’s.  22) 
to  solve  for  the  unknown  displacements.  The  quantities  [A]  and  [B]  or 
equivalently,  Ar>  and  Nlrs  are  calculated  once  based  on  the  equilibrium 
displacements.  Finally,  the  load  proportionality  parameter,  is 
incremented  until  a  sign  change  on  the  left  side  of  Eqn.  23  occurs, 
signifying  bifurcation. 


MATERIAL  DEGRADATION  INFLUENCES  ON  COMPOSITES 


Extreme  temperature,  near  or  beyond  the  polymer’s  glass  transition 
temperature,  and  moisture  (water  absorbed  by  the  polymer  resin)  can  cause 
severe  degradation  of  the  superior  properties  of  the  composites.  The 
degradation  factors  which  influence  the  resin  properties  are  increased 
with  temperature,  with  the  absorption  by  the  polymeric  resin  material  of 
a  swelling  agent  such  as  water  vapor,  and  with  the  sudden  expansion  of 
absorbed  gases  in  the  resin  [21].  The  composite  fibers,  on  the  other 
hand,  which  are  typically  graphite,  are  not  affected  by  either  moisture 
or  temperature.  Therefore,  the  resin  which  is  influenced  by  temperature 
and  moisture,  plays  the  most,  important  part  in  the  changing  of  the 
composite’s  properties. 

Temperature  and  moisture  act  on  the  polymeric  resin  in  two  ways. 

The  first  is  that  moisture  absorption  makes  the  resin  swell,  causing  a 
change  in  the  residual  stresses  of  the  composite  and  possibly  micro-crack 
formation.  This  matrix  swelling  and  the  rapid  heating  eventually  leads 
to  surface  cracks.  Secondly,  moisture  absorbed  in  the  resin  results  in 
plasticization  of  the  resin,  the  result  of  the  lowering  of  the  glass 
transition  temperature ,  T* .  The  glass  transition  temperature  is  actually 
a  temperature  range  below  which  the  resin  is  essentially  brittle  and 
above  which  it  behaves  rubbery.  As  the  resin  changes,  one  observes  a 
decrease  in  the  tensile  properties  and  a  reduction  in  the  shear  moduli  of 
the  composite  as  well  as  a  slight  increase  in  the  longitudinal  elastic 
moduli  [23], 


The  primary  mechanism  for  absorption  of  moisture  is  through 
diffusion.  In  this  analysis,  Pick’s  second  law  of  diffusion  [28],  which 
under  certain  circumstances  has  been  shown  to  correlate  with  test  data 
[23],  was  used.  Pick  developed  this  equation  in  1855  by  drawing  an 
analogy  between  heat  conduction  in  a  solid  and  diffusion  through  the 
solid.  Pick’s  equation  can  be  written  as 

ac/  at  =  k  (  a 2 c/  az1)  (25) 

where  C  is  a  measure  of  moisture  of  the  laminate  as  a  function  of  time 
and  distance  through  the  thickness,  z  equals  the  space  coordinate 
measured  normal  to  the  surface,  K  the  diffusion  constant,  and  t.  the  time 
C  is  expressed  as  a  ratio  of  the  gain  in  the  weight  of  the  laminate  due 
to  the  absorption  of  moisture  divided  by  the  original  weight  of  the 
laminate . 

The  solution  to  this  partial  differential  equation  (Eqn.  25)  with 
boundary  and  initial  conditions  pertinent  to  the  problem  is  shown  below. 
This  series  solution  in  a  slightly  different  form  is  found  in  Sec.  4.3.3 
of  Ref . [28] 

C ( z , t )  =  Ci  +  ( C 2  -  C i ) * ( z / h )  +  (26) 

oo 

(z/rr  )  ( ( C  2  cos  n  rr  -  Ci)/n)  +  (sin  ( n  7r  z/h  )  )  *  ( e  xp  ( -  K  n 2  rr2t)/h2)  + 

n  =  l 
30 

( 4Co  /  7r  )  ^  (l/2m+l)*(sin  ((2m  +  l)  irz)/h)*(exp  ( -K  ( 2m+ 1 ) 2  tt  2 1 ) /h2 ) 

m=0 


where  C  is  as  previously  defined,  Co  the  initial,  uniform  moisture 
concentration  through  the  thickness  of  the  laminate,  Ci  and  C2  the 
moisture  concentration  boundary  conditions  at  the  inside  (-2)  and  the 
outside  ( +  z )  surfaces  of  the  laminate,  respectively,  and  h  the  thickness 
of  the  laminate. 

This  series  solution,  which  is  a  combination  of  a  steady  state 
moisture  effect  and  a  transient  moisture  distribution,  allows  for 
the  determination  of  the  through  the  thickness  moisture  distribution. 
Because  of  a  series  solution,  accuracy  is  dependent  on  the  number  of 
terms  that  are  used.  Snead  [19]  wrote  a  computer  program  which 
calculates  the  solution  to  the  series  approximation  using  the  first  14 
terms  of  the  approximation,  and  thus  the  mechanical  properties  of  each 
ply  can  be  determined. 

There  are  limitations  on  the  application  of  Fick’s  equation.  The 
series  solution  of  Fick’s  equation  assumes  that  the  moisture  diffusion 
coefficient  K  is  constant  (for  this  study  K  equals  0.52537  x  10'10 
in2/sec  [24]).  In  reality,  the  diffusion  coefficient  is  a  function  of 
the  laminates  temperature  and  moisture  concentration  and  a  series 
solution  to  the  Fick  equation  using  this  functional  approach  can  be  found 
in  Ref.  [29].  However,  for  composite  material  studies,  the  diffusion 
process  may  be  assumed  to  take  place  at  a  constant  temperature.  This  is 
true  because  moisture  diffusion  is  a  relatively  slow  process  with  many 
months  and  years  required  before  the  moisture  concentration  distribution 
through  the  laminate  achieves  equilibrium.  Thus,  K  is  considered  to  be  a 
constant  as  is  usually  done  when  the  Fick  equation  is  incorporated  into 
composite  material  studies  [19,24], 


Fick’s  equation  to  model  moisture  diffusion  in  composites  is  also 
affected  by  rapid  temperature  changes.  Rapid  thermal  heating  of  the 
laminate  near  the  materials’  T*  has  been  found  to  increase  the  rate  of 
moisture  weight  gain  above  that  predicted  by  Fick’s  equation  [23].  This 
increase  is  believed  to  be  due  to  the  development  of  surface  crazing  arid 
cracking  brought  about  by  the  rapid  heating  arid  resin  swelling.  With  the 
restrictions  of  no  rapid  heating,  no  surface  crazing  or  cracking,  and 
assuming  K  to  be  constant,  Fick’s  equation  gives  a  good  initial 
approximation  of  the  moisture  concentration  distribution  through  a 
composite  laminate  [24]. 


III.  FINITE  ELEMENT  MODELLING 

Panel  Properties 

For  the  finite  element  analysis,  a  graphite  epoxy  composite  panel 
that  is  12  inches  long,  a  c ircumferen tial  length  of  12  inches,  and  a 
radius  of  curvature  of  12  inches  was  modeled.  The  panel’s  dimensions  and 
material  properties  are  shown  in  Fig.  3.1.  The  panel  consists  of  eight 
0.005  inch  plies  for  a  total  thickness  of  0.04  inches.  Two  different 
symmetric  ply  orientations  were  examined,  (0/+45/-45/90/90/-45/+45/0) 
referred  to  as  (0/+45/-45/90)s  and  (0/-45/+45/90)s .  All  the  panels  were 
loaded  along  the  top  edge  which  was  clamped  with  the  u  displacement  free. 
The  vertical  edges  of  the  panel  had  u  and  v  displacements  (x  and  y 
direction  displacements,  respectively)  and  rotation  about  the  x  axis,  Ru , 
free.  The  bottom  of  the  panel  was  clamped  with  no  free  degrees  of 


freedom  (DOF). 


ss 
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DIMENSIONS  AND  MATERIAL  PROPERTIES 

MATERIAL:  GRAPHITE-EPOXY 
RADIUS:  R-12” 

LENGTH:  L-12” 

NUMBER  OF  PLIES:  8 

ORIENTATION  OF  PLIES, 0:  (0/45/-45/90)* 

(0/-45/+45/90)s 

THICKNESS.  8  PLIES  AT  0.005”  *  0.04" 

ELASTIC  MODULI:  E,  -18850  ksi,  E--1413.8  ksi 

SHEAR  MODULUS:  G-855  ksi 

POISSON'S  RATIO:  v, 2 -0.3 

x.  y,  v  STRUCTURAL  COORDINATES 

u,  v,  w:  DISPLACEMENTS 

Ry,  \  R^  ROTATIONS 


Fig.  3.1  Panel  Notation  and  Material  Properties 


Element  Selection 


For  the  finite  element  analysis,  the  element  selected  is  a 
quadrilateral  plate  element  referred  to  as  the  9UAF  411  in  the  STAGSC-1 
documentation  [26].  The  QUAF  411  consists  of  three  translational 
displacements,  two  in-plane  rotations  and  two  independent  normal 
rotations  at  each  corner  node,  plus  tangential  displacements  at  each  of 
the  four  midside  nodes  for  a  total  of  32  DOF.  QUAF  411,  shown  in  Fig. 
3.2,  was  developed  to  remove  the  displacement  incompatibilities  which 
exist  when  flat  shell  elements  do  not  lie  in  the  same  plane.  The 


in-plane  shape  functions  are  represented  by  cubic  polynomials  and  if  the 
edge  displacements  are  to  be  compatible,  the  transverse  displacement 
shape  functions  must  be  of  the  same  order.  This  is  accomplished  by  the 


use  of  cwo  normal  rotations  at  each  corner  node  and  tangential 
displacements  at  mid-side  nodes.  The  difference  between  these  two 
rotations  yields  a  shear  strain  at  each  corner  node  which  is  introduced 
as  an  extra  degree  of  freedom.  The  shape  functions  for  the  u  and  v 
displacements  within  the  QUAF  411  element  are  cubic  polynomials  parallel 
to  the  edge  and  quadratic  perpendicular  to  the  edge.  Also,  the  bending 
shape  functions  are  cubic  in  both  directions. 

This  thesis  requires  a  nonlinear  collapse  analysis  to  be  carried  out 
by  the  STAGSC-1  computer  program  and  the  QUAF  411  element  was  the  best 


available  element  to  use  [26].  The  drawback  to  the  QUAF  411  element  is 
that  there  is  no  transition  element  in  the  STAGSC-1  library  so  grid 
selection  becomes  very  important. 
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Grid  Sizing 


Grid  selection  is  a  trade-off  between  cost  and  accuracy.  Lee  [18] 
pointed  out  that  the  elements  should  be  as  close  to  an  aspect  ratio 
(axial  dimension  divided  by  circumferential  dimension)  of  unity  as  can 
possibly  be  achieved.  Elements  with  an  aspect  ratio  of  0.5  should  be  no 
closer  than  at  least  two  inches  from  the  cutout  and  elements  with  an 
aspect  ratio  of  2  should  be  no  closer  than  at  least  one  and  one-half 
inches  from  the  cutout  considering  cutout  areas  of  up  to  112  of  the 
panels  surface.  These  important  parameters  were  used  and  an  element  mesh 
was  constructed.  The  mesh  used  in  this  paper  was  refined  around  the 
cutout  to  obtain  better  accuracy. 

The  refined  mesh  that  was  chosen  for  the  analysis  was  compared  to 
the  17  x  19  grid  of  Lee  [13]  for  a  centrally  located  2"  x  2"  cutout.  The 
major  difference  between  the  two  mesh  models  is  located  in  the  vicinity 
of  the  cutouts  reentrant  corners.  The  mesh  used  herein  is  302  finer  in 
these  areas.  Shown  in  Table  3.1  are  the  results  of  the  comparison  of  the 
17  x  19  grid  of  Lee  with  the  refined  17  x  20  grid  used  in  this  analysis. 
As  the  table  shows,  the  17  x  20  grid  gives  excellent,  results.  Therefore, 
since  the  grid  refinement  chosen  compares  favorably,  it  was  used  for  the 
cutout  and  material  degradation  analysis. 


3.1:  Mesh  Refinement  Results  for  a  2"  x  2"  Cutout 
(Note-  a  different  set  of  material  properties 
was  used  to  match  the  published  results) 

Published  results  [18]  Refined  Grid 

17  x  19  grid  17  x  20  grid 

2381  Active  DOF  255b  Active  DOF 

08.1  lbs/in  collapse  load  208.7  lbs/in  collapse  load 
26  incremental  steps  24  incremental  steps 


0.3Z  difference 


IV.  DISCUSSION  AND  RESULTS 


The  main  study  area  of  this  thesis  is  the  material  degradation 
effects  considering  cutout  location.  Thus,  in  order  to  establish  the 
critical  characteristics  of  the  effect  that  a  true  geometric  imperfection 
(cutout)  has  on  a  composite  panel,  a  study  of  the  effect  of  the 
discontinuity  location  in  the  panel  has  been  carried  out.  The  cutouts 
consisted  of  a  2"  x  2"  cutout  which  removed  2.732  of  the  panels  total 
surface  area  and  the  4"  x  4”  cutout  which  removed  11.12  of  the  total 
area.  The  different  ply  layups,  different  size  cutouts  and  the  cutout 
locations  studied  are  shown  in  Table  4.1.  The  reader  must  recognize  that 
the  results  presented  within  this  thesis  are  based  on  a  finite  element 
study  and  no  attempt  has  been  made  to  compare  the  results  with 
experimental  values. 

The  first  cutout,  considered  was  located  in  the  middle  of  the  panel, 
cutout  location  A  (Fig.  1.2).  A  collapse  analysis  was  completed  on  two 
symmetric  ply  orientations,  (0/+45/-45/90)s  and  (0/-45/+45/90)s  (collapse 
is  considered  the  point  in  which  the  load  can  no  longer  increase  with 
increased  displacement).  The  different  orientations  had  no  effect  on  the 
collapse  load  of  199.4  lbs/in.  However,  there  was  a  difference  in  the 
radial  displacement  contours  as  shown  in  Figs.  4.1  arid  4.2.  As  can  be 
seen  in  Fig.  4.1,  at  the  final  load  level  (collapse  load),  one  can  see 
the  loss  of  symmetry  and  the  formation  of  a  trough  of  radial  displacement 
that  runs  from  the  upper  left  of  the  panel  to  the  lower  right.  In 
looking  at  Fig.  4.2,  (for  the  (0/-45/+45/90)s  orientation)  the  trough  of 
radial  displacement  runs  from  the  upper  right 
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to  the  panel’s  lower  left.  This  difference  in  the  orientation  of  the 
troughs  seems  to  be  related  to  the  orientation  of  the  second  ply  (taken 
from  the  outer  surface  edge).  A  (0/+45/-45/90) s  orientation  produces  a 
trough  along  the  direction  of  the  +45°  ply  (ply  angle  defined  in  Fig. 

3.1)  while  a  (0/-45/r45/90)s  orientation  produces  a  trough  along  the 
order  of  the  -45°  ply. 

In  order  to  explain  the  unexpected  collapse  displacement,  for  the 
2"  x  2"  cutout,  a  comparison  was  made  between  the  eigenvector  solution 
and  the  radial  contour  plot.  Fig.  4.3  represents  the  eigenvector  for  the 
( 0 /  +  4 5 / -45/90) s  orientation,  and  one  should  observe  that  a  trough  similar 
to  the  panel’s  collapse  mode  is  appearing.  However,  the  overall 
eigenvector  shape  does  not  resemble  that  of  the  collapse  mode.  Further 
study  into  this  area  was  required.  Fig.  4.3.1  represents  the  total 
linear  solution  prior  to  the  bifurcation  (buckling)  point  and  Fig.  4.4 
represents  the  nonlinear  2“  x  2"  collapse  analysis  at  502  of  the  total 
collapse  force.  As  the  two  plots  show,  the  displacement  contours  are 
nearly  identical.  This  clearly  shows  that  the  2"  x  2"  nonlinear  analysis 
closely  followed  the  linear  solution  through  the  initial  stages  of 
loading.  Further  proof  of  this  linear  phenomena  can  be  found  in  Fig. 

4.5,  the  plot  of  the  load  displacement  curves.  As  the  plot  shows,  the  2" 
x  2"  cutout  closely  followed  the  linear  bifurcation  plot..  There  is, 
however,  a  difference  between  the  final  bifurcation  load  point  and  the 
collapse  load  of  roughly  132.  As  the  linear  prebuckling  contour 
indicates,  a  drastic  change  in  the  displacement  function  occurs  at 
bifurcation  in  which  the  trough  is  developed.  This  shows  that  a  greater 
force  is  needed  to  go  from  the  linear  displacement  contour  plot  to  the 
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g.  4.4  Radial  Displacement  Contour  Plot  for  2"  x  2'  Cut 
of  Collapse  Load,  101  lbs/in,  (0/+45/-45/90)s ,  C 
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eigenvector  then  the  force  needed  to  collapse  the  panel,  as  would  be 


expected.  This  is  primarily  true  because  the  nonlinear  collapse  analysis 
contour  plot  responds  gradually  to  the  trough  development.  The  collapse 
analysis  is  still  a  nonlinear  phenomena  even  though  it  appears  to  follow 
the  linear  path.  Shown  in  Fig.  4.5.1  is  the  load  displacement,  curve  for 
a  point  located  halfway  along  the  left  side  of  the  cutout.  This  curve 
shows  that  the  solution  does  become  nonlinear  as  collapse  approaches.  It 


can  also  be  shown  using  this  same  analysis  that  for  the  (0/-45/+45/90) s 
panel  the  eigenvector  is  such  that,  as  mentioned  previously,  the 
displacement  trough  is  directed  in  an  orientation  taken  on  by  a  negative 


angle,  Fig.  4.6. 

A  further  description  can  be  made  relative  to  the  radial 
displacement  function.  As  pointed  out  previously,  there  was  an  initial 
surprise  in  seeing  unsymmetric  deflection  movements.  The  symmetry  was 
evident  when  the  panel  is  looked  at  for  lower  load  values  or  what  is 
called  the  linear  portion.  In  this  linear  portion,  the  different 


orientations,  (0/+45/-45/90)s  and  (0/-45/+45/90)s ,  did  not  matter  in  the 
solutions  displacement  contours  as  they  were  identical  and  symmetric.  At 
the  higher  loadings  the  influence  of  the  cutout  comes  into  play  and  the 
final  contour  plot  results  in  the  unsymmetrical  secondary  loading 
arrangement  (the  eigenvectors)  as  previously  mentioned.  The  question  may 
be  asked  as  to  what  feature  of  the  stiffness  coefficients  could  indicate 
an  unsymmetrical  and  different  displacement  function  within  the  shell  by 
only  changing  the  ply  orientation  (actually  just  switching  the  second  and 
third  plies)  of  the  quasi-isotropic  layup.  The  answer  lies  in  the  Di6 
and  Die  terms  of  the  bending  stiffness  parameter  matrix  since  the  other 


i 


stiffness  parameters  are  equal.  These  Du  coefficients  are  related  to 
the  twisting  moment,  Mx>.  The  Mxy  equation  is  written  as 

Hxy  -  Dlt  +  D2i  *y  t  Dit  *xf  (27) 

where  the  *’s  are  the  changes  in  the  surface  curvatures  and  the  Du’s  are 
the  bending  stiffness  elements  [24,25].  In  the  (0/  +  45/-45/90 ) s 
arrangement,  the  Du  and  Du  are  positive  and  in  the  (0/-45/  +  45/'’0)s  ply 
orientation  the  terms  have  the  same  magnitude  but.  are  negative  in  sign. 
This  means  that  Mxy  will  change  with  ply  orientation  and  since  this 
moment,  function  dominates  the  collapse  for  the  small  2"  x  2"  cutout, 
different  contour  arrangements  are  produced.  Thus,  there  is  a  tie-in  with 
the  twisting  moment,  the  curvature  and  the  displacement  w.  However,  in 
the  linear  loading  portion,  since  there  is  very  little  bending  moment  the 
Au  (extensional  stiffness  parameters)  matrix  dominates  the  solution  and 
Du  and  Du  can  not  affect  the  displacement  contours.  Therefore,  at  the 
lower  loads  (linear  range)  the  displacement,  contours  do  not  follow  the 
second  ply  angle.  In  looking  at  the  effect  a  4 "  x  4 "  cutout  had  on  the 
panel,  one  notices  a  difference  from  the  2"  x  2”  cutout.  First,  it  can 
be  noticed  in  Fig.  4.5,  collapse  load  versus  end  shortening  and  in  Fig. 
4.5.2,  load  displacement,  curve  of  a  point,  located  halfway  along  the  left 
side  of  the  cutout,  that  the  solution  is  totally  nonlinear  for  the  4"  x 
4"  cutout.  This,  as  previously  mentioned,  is  riot  true  when  considering 
the  2"  x  2"  cutout.  Furthermore,  it  can  be  observed  that  collapse  load 
and  bifurcation  load  are  very  close.  It  becomes  obvious  that  the 
non  1  inear l ty  surrounding  the  4"  x  4"  cutout  is  associated  with  a  larger 


amount  of  bending  which  allows  the  collapse  to  seek  a  level  similar  to 
the  linear  load.  Yet.,  this  same  phenomena  is  not  true  for  the  2"  x  2" 
cutout  since  it  is  depicting  a  linear  relationship  in  the  area 
surrounding  the  actual  cutout.  Linearity  is  not  conducive  to  direct 
bending  moment,  such  as  Mx  and  My .  Thus,  the  separation  between 
bifurcation  and  collapse  load  is  apparent  when  one  considers  the  2"  x  2" 
cutout. 

Once  again,  looking  at  the  larger  cutout  (4“  x  4"),  the  radial 
displacement  contours  are  investigated.  Fig.  4.7  represents  the  contour 
displacements  for  a  collapse  load  of  90.0  lbs/iri  and  a  ply  orientation  of 
( 0 /  +  4 5 / - 4 5 / 9 0 )  s  .  For  a  ply  orientation  of  (0/-45/  +  45/90)s  (plot  not 
shown)  the  same  symmetric  contour  displacements  can  be  seen  with  only  a 
slight  amount  of  non-symmetry  that  can  be  observed  at  the  center  of  the 
cutout’s  vertical  edges.  In  looking  at  the  contour  plot  one  can  see  that 
the  displacement  pattern  at  collapse  is  more  symmetric  about  the 
circumferential  and  longitudinal  axes  than  the  similar  panel  with  a 
smaller  cutout.  This  is  because  the  larger  cutout  produces  the  direct 
bending  effect  earlier  and  the  moment  change  is  more  spread  out, 
therefore,  a  more  symmetric  displacement  pattern.  It  can  also  be  seen 
that  the  larger  (4"  x  4”)  cutout  does  not  display  the  inward  displacement 
trough  along  the  second  ply  angle  as  was  seen  in  the  smaller  (2"  x  2") 
cutout  giving  evidence  the  Mxy  function  is  not  as  dominant  relative  to 
the  other  moment  resultants. 

Fig.  4.8  is  the  radial  component  eigenvector  contour  plot  of  the 
( 0/ +45/ -45/90 ) s  panel  with  the  larger  cutout  (4"  x  4").  In  comparing  the 
eigenvector  plot  to  the  displacement  contours  at  the  collapse  load,  there 
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Fig.  4.8 


Eigenvector  Contour  Plot  for  4' 
(0/*45/-45/90)s ,  Cutout  4 


x  4*  Cutout 
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is  very  little  similarity.  The  radial  displacement  contour  plot.  Fig. 
4.7,  at  collapse  produces  a  symmetric  displacement  along  a  vertical 
centerline.  On  the  other  hand,  the  eigenvector  plot,  Fig.  4.3,  has  a 
displacement  pattern  that  is  slightly  tilted  and  the  overall  pattern  does 
not  resemble  the  collapse  displacement  pattern.  Fig.  4.3.1  represents 
the  linear  solution  prior  to  bifurcation  and  closely  resembles  the 
nonlinear  collapse  analysis  plot,  Fig.  4.7.  In  essence,  when  the  panel 
collapsed  it  did  not  go  into  the  shape  of  the  eigenvector  as  the  panels 
with  the  smaller  cutouts  did.  If  the  panel  absorbs  energy  primarily 
through  its  axial  stiffness  (smaller  cutout  and  uncut  panel)  then  the 
displacement  field  at  collapse  is  eigenvector  oriented.  However,  if  the 
panel  absorbs  energy  through  both  the  axial  and  bending  stiffness  then 
the  displacement  field  is  not  eigenvector  oriented.  This  difference  can 
also  be  explained  by  comparing  the  amounts  of  radial  displacement. 

In  looking  at  the  radial  displacements,  the  small  cutouts  had 
maximum  radial  displacements  approximately  equal  to  the  panel  thickness 
(0.04  inches).  The  largest  radial  displacement  in  the  panels  with  small 
cutouts  was  0.047  inches  and  that  was  on  the  panel  of  ( 0 /  +  4  5 / -45/30 )  s .  On 
the  other  hand,  the  panels  with  the  larger  cutout  (4"  x  4")  had  radial 
displacements  of  0.0t.4  inches  which  is  over  one  and  one-half  the  panels 
thickness.  Therefore,  the  displacement  field  generated  by  the  2“  x  2“ 
cutout  is  not  going  to  have  that  great  of  an  effect  on  the  nonlinear 
collapse  mode  since  the  radial  displacements  are  small  (near  the  panel’s 
thickness).  Thus,  at  collapse  the  panel  will  behave  almost  linear  and 
will  have  the  shape  of  the  linear  bifurcation.  Since  the  radial 
displacements  for  the  panel  with  the  large  cutout  are  one  and  one-half 
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Fig.  4 . 


.1  Linear  Radial  Displacenent  Contour  Plot  for  4'  x  4" 
Cutout  Prior  to  Bifurcation  Buckling,  (0/+45/-45/90) 
Cutout  Position  A. 


the  panels  thickness,  their  effect  on  the  nonlinear  collapse  analysis 
will  be  greater.  With  these  large  displacements  (greater  than  the 
panel’s  thickness)  the  nonlinear  effect  of  the  cutout,  on  the  collapse 
analysis  dictates  the  collapse  displacement  pattern  as  discussed 
previously.  This  pattern  is  symmetric  and  there  is  no  shift  to  the 
secondary  loading  path.  The  nonlinearity  in  the  presence  of  the  4"  x  4'' 
cutout  is  felt  by  the  panel  throughout  its  loading  history  while  the 
nonlinearity  of  the  smaller  2"  x  2“  cutout  does  not  effect  the  panel’s 
displacement  pattern  as  much. 

To  continue  on  with  the  cutout  placement  study,  a  ( 0/+45/-45/90) s 
ply  orientation  2''  x  2"  cutout  was  placed  near  the  top  of  the  panel 
(cutout  location  B).  This  cutout  location  produced  a  collapse  load  of 
201.15  lbs/in,  or  a  0.3?  change  from  the  2"  x  2"  cutout  located  at  the 
exact  middle  (cutout  location  A).  The  radial  displacement  contour  plot, 
Fig.  4.9,  looks  exactly  like  the  middle  cutout  (location  A)  except  that 
the  cutout  is  displaced  upward.  Fig. 4. 10  is  the  eigenvector  plot  for  the 
cutout  located  on  the  top  of  the  panel  (location  B).  Once  again  we  see 
that  at  collapse  the  contour  plot  has  shifted  into  the  secondary  loading 
path  ! the  eigenvector).  Since  the  top  cutout  acts  very  nearly  the  same 
as  the  middle  cutout,  its  analysis  would  yield  the  same  interpretation  as 
the  middle  c  u  L  ^ ...  t  1 1  s  c  u  s  s  e  d  earlier’. 

Next,  if  the  2"  x  2"  cutout  is  placed  along  one  side,  such  as  at 
location  C,  there  is  a  unique  dependence  on  the  ply  orientation.  A  ply 
orientation  of  (0/+45/-45/90)s  produces  a  collapse  load  of  147.18  lbs/in 
or  a  26.2?  decrease  from  the  2"  x  2"  middle  cutout.  With  a  ply 
orientation  of  ( 0/ -45/ +  45/90 ) s ,  a  collapse  load  of  178.6  lbs/  in  is  found, 
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g.  4.9  Radial  Displaceaent  Contour  Plot  for  2*  i  2*  Cutout  at 
Collapse  Load  of  201.2  lbs/in  (0/+45/-45/90)s ,  Cutout  8 
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Fig.  4.10 


Eigenvector  Contour  Plot  for  2’ 
(0/+45/-45/90)* ,  Cutout  B 


2"  Cutout 


which  is  10.42  less  than  the  middle  cutout  collapse  load  arid  17.62  more 
than  the  ( 0 /  +  4 5 / - 4 5 / 9 0 ) s  side  cutout.  Radial  displacement  contour  plots 
are  different  for  the  ( 0/  +  45/-45/90) s  and  the  ( 0/-45/  +  45/90 ) s ,  Fig.  4.11 
and  4.12  respectively.  In  looking  at  Fig.  4.11,  the  radial  displacement 
contour  plot  for  (0/ +45/-45/90) s ,  one  observes  again  the  inward 
displacement  trough  along  a  line  near  the  second  ply  angle  (+45°).  This 
trough  runs  along  the  same  angle  as  that  of  the  middle  cutout  (cutout  A). 
A  big  difference  is  noticed  when  Fig.  4.12  is  examined.  This  contour 
plot,  if  it  were  similar  to  the  middle  (0/-45/+45/90) s ,  would  have  a 
trough  along  the  -45°  degree  angle.  However,  there  is  no  trough  along 

the  ~ 4 5 3  degree  angle  as  we  had  seen  earlier.  What  is  seen  is  a 

semicircular  shaped  trough  that  starts  in  the  upper  left  and  ends  in  the 
lower  left.  The  semicircle  trough  allows  for  a  greater  curvature  and  the 
panel  absorbs  more  energy  and  achieves  a  higher  collapse  load.  Figs. 

4.13  and  4.14  represent  radial  displacement  eigenvector  contour  plots  for 
the  ( 0/ +45/  —  45/90 ) s  and  (0/-45/  +  45/90)s  respectively.  In  reviewing  Fig. 
4.11,  one  can  see  that  at  collapse  the  (0/+45/-45/90) s  ply  orientation 
has  shifted  to  the  secondary  loading  path  (the  eigenvector  shown  in  Fig. 
4.13).  However,  in  reviewing  Fig.  4.12  and  4.13  both  plots  of  the 
(0/-45/+45/90)s,  one  sees  that  at  collapse  the  displacement  contours  are 
not  in  the  secondary  (eigenvector)  mode.  This  associated  phenomena  can 
be  explained  again  by  looking  at  Mxy.  A  side  cutout  means  that  there  is 

less  material  between  the  cutout  and  the  edge.  Therefore,  loads  build  up 

at  a  higher  rate  between  the  cutout  and  the  side  edge.  Since  the  cutout 
is  closer  to  the  edge,  the  values  of  Dit  and  D26  (Eqn.  27)  and  there 
positive  or  negative  values  depending  on  ply  orientation,  determine  the 
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Fig.  4.11  Radial  Displacement  Contour  Plot  for  2*  i  2’  Cutout  at 
Collapse  Load  of  147.2  lbs/in  (0/*45/-45/90)s ,  Cutout  C 
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Fig.  4.12  Radial  Displacement  Contour  Plot  for  2*  *  2*  Cutout  at 
Collapse  Load  of  178.6  Ibs/in  (0/ -45/ +45/90)* ,  Cutout  C 
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Fig.  4. 


14  Eigenvector  Contour  Plot  for  2* 
(0/-45/  +  45/90)* ,  Cutout  C 


2*  Cutout 


amount  of  twisting  moment,  Mxy  on  the  panel.  The  ( 0/  +  45/ -45/^0 )  s 
orientation,  as  seen  earlier,  has  positive  values  for  Die  and  D 2 1  and 
thus  a  greater  twisting  moment  dominance  and  thus  a  lower  collapse  load. 
The  opposite  is  true  for  the  ( 0/-45/ +45/^0) s  orientation.  Here,  since 
the  D 1  &  and  D  2  6  values  are  negative,  the  twisting  moment  is  less  and  the 
panel  can  withstand  a  higher  collapse  load  then  the  (0/  +  45/-45/"^0)  s 
configuration. 

Load  displacement  plots  for  cutout  A  and  cutout  C  are  shown  in  Fig. 
4.15.  This  plot,  demonstrates  the  displacement  in  the  direction  of  the 
lead.  The  slight  curvature  near  the  collapse  load  expresses  the  fact 
that  the  problem  is  nonlinear.  Since  cutout  C  has  a  greater  slope  then 
cutout  A,  cutout  C  has  greater  axial  stiffness. 

Now  a  4"  x  4"  cutout  was  placed  along  one  side,  cutout  location  C. 
The  4"  x  4"  cutout  does  not  have  a  unique  dependence  on  the  ply 
orientation  as  was  seen  with  the  2"  x  2"  cutout.  A  ply  orientation  of 
( 0  /  +  45/ -45/90) s  produces  a  collapse  load  of  66.0  lbs/  in  or  a  2b.  11 
decrease  from  the  4"  x  4"  middle  cutout.  Now,  with  a  ply  orientation  of 
( 0/ - 4 5 /  +  4 5/ 9  0 ) s ,  a  collapse  load  of  66.0  lbs/ in  was  also  obtained.  The 
radial  displacement  contour  plots  are  the  same  for  the  (0/+45/-45/w0)s 
and  the  ( 0/-45/  +  45/90) s  ,  Fig.  4.16  and  4.17  respectively.  In  looking  at 
Fig.  4.16,  the  radial  displacement  contour  plot  for  ( 0/ M5/ -45/,J0 )  s ,  one 
observes  that  the  inward  displacement  trough  is  not  along  the  second  ply 
angle  (+45°).  The  same  displacement  are  noticed  when  Fig.  4.17  is  also 
examined.  As  was  the  case  with  the  middle  4"  x  4"  cutout,  there  is  more 
symmetry  about  the  circumferential  and  longitudinal  axes  then  similar 
panels  with  cutouts.  As  mentioned  previously,  the  larger  cutout  produce 
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the  direct  bending  effect  earlier  and  the  moment  change  is  more  spread 
out,  therefore,  a  more  symmetric  displacement  pattern.  As  before,  when 
the  panels  collapsed,  they  did  not  go  into  the  shape  of  the  eigenvector 
(plot  not  shown)  as  the  panels  with  the  smaller  cutouts  did. 

This  same  phenomena  that  was  observed  for  the  cutout  at  location  C 
can  be  seen  for  a  cutout  located  along  the  right  side,  cutout  location  D 
(Fig.  1.2).  There,  however,  is  one  slight  change.  Fig.  4. IS  the  cutout 
D  ( 0/ r45/ -45/90 ) s  ply  orientation  contour  plot  exhibits  exactly  the  same 
characteristics  as  the  (0/-45/  +  45/90)s  ply  orientation  plot.,  Fig.  4.12, 
for  cutout  location  C.  The  same  is  true  for  the  cutout  D  ( 0/ -45/ +45/^0 ) s 
plot,  Fig.  4. 1'5,  and  the  cutout  C  (0/  +  45/-45/90)s  plot,  Fig.  4.11.  That 
is,  there  is  a  swapping  phenomena.  The  +45°  second  ply  arrangement  on 
the  right  side  looks  like  and  has  the  same  collapse  load  as  that  of  the 
-45°  second  ply  arrangement  of  the  left  side  cutout.  Likewise,  the  -45° 
right  side  arrangement  acts  exactly  like  the  +45°  left  side  arrangement. 
Due  to  this  special  phenomena,  the  analysis  of  the  right  side  would  be 
exactly  the  same  as  its  counterpart,  the  left  side,  as  discussed  earlier. 

The  results  of  the  discontinuity  location  study  is  presented  in 
Table  4.2.  For  the  small  cutout  (2"  x  2"),  the  different  ply 
orientations  really  do  not  matter  for  cutouts  located  on  the  vertical 
axis  of  symmetry.  Flowever,  if  a  cutout  is  repositioned  to  either  side  of 
this  vertical  axis  of  symmetry  the  ply  orientation  makes  a  significant 
difference.  For  the  larger  cutout  (4”  x  4")  the  ply  orientation  does  not 
make  a  difference  as  the  problem  has  become  totally  nonlinear. 
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Fig.  4.18  Radial  Displacement  Contour  Plot  for  2*  i  2' 
Collapse  Load  of  178.6  Ibs/in  (0/+45/_45/90)s 
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Table  4.2:  Discontinuity  Location  Study  Results 


Cutout 


Collapse 


Linear 


percent  of 


Location  Orientation  Load  (lbs/in)  Bifurcation  (lbs/in)  Collapse  at.  A 


( 0/ +45/ -45/90 ) 5 

199.4 

229.014 

10/ -45/ +45/90) s 

199.4 

229.014 

( 0/  +  45/ -45/90 ) > 

201.15 

222.035 

(  0/ -45/ +45/90  )  s 

201.15 

222.035 

( 0/  +  45/  —  45/90 ) > 

147.18 

198.479 

(0/  —  45/ +45/90) 3 

178.6 

193.530 

( 0/  +  45/ -45/90) s 

173.6 

198.530 

( 0/ -45/ +45/90 )  3 

147.16 

193.479 

4 "  x  4 " 


(  0/  +  45/-45/90) * 
(0/-45/M5/*0). 


■=>5.32 

95.32 

32.57 


(0/+45/-45/90) 
( 0 / - 4  5 / +45/90) 


66 . 0 


32.57 


Material  Degradation  Effects  Study 


The  previous  study  indicated  that  the  critical  cutout  locations 
for  establishing  upper  and  lower  bound  effects  can  be  attributed  to  the 
center  and  eccentric  positions.  Thus,  these  two  positions  were 
investigated  relative  to  material  degradation  effects.  For  the  material 
degradation  investigation  carried  out  within  this  study,  a  combination  of 
moisture  conditions,  temperature  and  time  intervals  were  looked  at  in 
order  to  produce  the  widest  change  in  phenomena  possible.  The  weight 
gain  due  to  moisture  concentration  distribution  and  the  moisture  and 
temperature  induced  degradations  in  the  E2  and  G12  moduli  are  shown  in 
Fig.  4.20  and  4.21,  respectively  [24].  For  this  analysis  two  temperatures 
were  examined,  80  °F  (300  °K)  and  250  °F  (394  °K).  The  80  °F  temperature 
represents  room  temperature  conditions  and  the  250  °F  temperature  is 
roughly  100  degrees  below  the  glass  transition  temperature.  This  250  °F 
temperature  does  represent  the  point  where  the  material  properties  really 
start  to  play  an  important  role.  Two  moisture  concentration  conditions 
were  chosen  (Table  4.3).  One  was  symmetric  (C«-0.0,  Ci =0.0105, 

C2  -0 . 0105 )  and  the  other  urisymmetric  (Co=0.0,  C 1  =0.0,  02  =  0.0105).  The 
moisture  concentration  is  measured  as  a  percentage  of  the  original  weight 
gained  through  moisture  absorption.  For  the  AS4/3501  graphite/epoxy,  the 
saturation  moisture  concentration  is  assumed  to  be  1.05Z.  Finally,  two 
different  time  events  were  chosen,  a  short  0.35  days  ( t.  *  =  0 . 0  0 1 )  and  a 
steady  state  176.24  days  ( t *  =  0 . 5 )  where  t*  is  a  normalized  time  equal  to 
K t / h 2 .  Shown  in  Table  4.4  are  the  relations  between  real  and 


dimensionless  time.  With  the  combinations  mentioned  above,  there  will  be 
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WEIGHT 


DISTANCE  THROUGH  THICKNESS  (in.) 


Fig.  4.20  Moisture  Concentration  Distribution  for  Unsytaetric  and 
Syaaetnc  Moisture  Conditions 
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Fig.  4.21  E*  and  Gu  Degradation  vs  Teeperature  at  Constant 
Values  of  Moisture  Concentration 
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Table  4.3 


»! 


Moisture  Conditions 

Moisture  Conditions 


Con d  .  No  • 

CO 

Cl 

c2 

1 

0.00 

0.00 

0.0103 

2 

0.00 

0.01  OS 

0.0105 

Table  4.4  Relation  Between  Real  and  Disensionless  Tiae 


Relation 

Between  Real  and 

Dimensionless  Tlae 

Real  Time 

Real  Time 

Diaensionless  Time 

(aec) 

(days) 

0.0 

0.0 

0.0 

3.045E04 

0.35 

0.001 

3.04  5E05 

3.  52 

0.01 

3.04  5E0t 

35.2- 

0.  1 

1.527EG7 

1 7fc .  24 

0.  5 

These  Claes 

were  calculated 

using  K  »  0. 52537E- 

.V--.  •-  ■-  -- 


a  total  of  3  runs  for  each  case.  Two  cases  were  chosen  for  this  study,  a 
2"  x  2"  cutout  located  in  the  middle  of  the  panel  (cutout  location  A)  and 
a  2"  x  2"  cutout  located  on  the  side  of  the  panel  (cutout  location  C) 
which  were  found  to  represent  the  overall  nature  of  cutout  location. 

The  easiest  way  to  analyze  the  material  degradation  effects  is 
through  the  use  of  Figs.  4.22  and  4.23.  In  these  plots  Nxong  represents 
the  bifurcation  load  for  the  room  temperature  condition  at  time  equal  to 
zero  ( t*  =  0 . 0 )  for  the  various  cutout  conditions  investigated.  This 
condition  is  unaffected  by  either  temperature  or  moisture  degradations. 

As  was  expected,  the  panel  bifurcation  load  decreased  with  increasing 
temperature  and  absorbed  moisture.  Fig.  4.22,  the  unsymmetric  moisture 
condition  plot,  shows  that  as  time  increases  there  is  little  effect  on 
the  reduction  in  collapse  load  for  the  cutouts  studied.  Whereas,  for  the 
panel  with  no  cutout,  there  is  a  significant  reduction  in  the  collapse 
load  as  the  temperature  is  increased.  Moreover,  there  is  very  little 
change  in  collapse  for  the  two  cutout  locations  considered  as  time 
progresses.  In  Fig.  4.23,  the  symmetric  moisture  condition  plot,  there 
is  an  effect  on  the  reduction  in  the  collapse  load  as  temperature  is 
increased  for  both  the  non-cutout  and  cutout  panel.  Again,  there  is  very 
little  effect  due  to  different  cutout  locations  on  the  change  in  the 
reduction  of  the  collapse  load  due  to  the  material  degradation  effects. 

Looking  once  again  at  Figs.  4.22  and  4.23,  involving  buckling  versus 
time  for  two  different  boundary  conditions,  one  observes  certain 
characteristics  (previously  mentioned),  especially  the  fact  that  the 
symmetric  boundary  condition  produces  a  greater  reduction  of  the  buckling 
load  as  temperature  and  time  are  increased.  In  order  to  explain  these 


characteristics,  look  at  Figs.  4.20  and  4.21.  For  the  unsymmetr ic  case 


at  t*-0.001,  the  weight  gain  due  to  moisture  (Fig.  4.20)  conditions  is 
concentrated  on  only  one  edge  of  the  panel.  In  comparison,  the  symmetric 
moisture  conditions  produce  twice  as  much  weight  gain.  If  one  now  looks 
at  Fig.  4.21,  it  can  be  seen  that  the  moisture  concentration  percentages 
follow  a  negative  slope.  The  symmetric  moisture  condition,  which  is 
twice  as  concentrated  as  the  unsymmetnc  case,  would  then  produce  greater 
changes  in  the  material  properties  as  temperature  increases.  If  the  same 
comparison  is  made  for  t *  :<J .  5  ,  an  even  greater  change  can  be  seen. 
Therefore,  the  symmetric  boundary  condition,  which  is  affected  the  most, 
produces  a  greater  affect  on  the  reduction  of  the  buckling  load  as  time 
and  temperature  are  increased. 

The  radial  displacement  contour  plots  for  the  temperature  and 
moisture  affected  panels  produce  the  same  results  obtained  for  the 
non-material  degraded  panels.  There  is,  however,  a  reduction  in  the 
collapse  load  as  the  moisture  and  temperature  changes  take  effect.  Shown 
in  Figs.  4.24  and  4.25  are  the  radial  displacement  contour  plots  for  the 
2"  .<  2"  cutout  A  at  collapse  load  for  room  temperature  (SO  °F)  and  for 
the  high  temperature  (250  0 F i  both  at  symmetric  moisture  conditions 
respectively.  As  can  be  seen  by  both  Figs.  4.24  and  4.25  the  conditions 
at  collapse  are  almost  identical  for  both  cases  as  well  as  looking  like 


the  non-affected  (no  temperature  or  moisture  degradation)  panels.  Fig. 
4.26  represents  the  material  degradation  effect  of  the  4"  x  4"  cutout  A 
which  has  nearly  the  same  contour  displacement  plot  as  the  normal 
non-affected  4"  x  4"  panel. 


. safc 


Now  shown  in  Table  4.5  are  the  results  of  the  material  degradation 
study.  This  table  clearly  shows  that  the  symmetric  moisture  condition 
produces  a  greater  affect  on  the  collapse  load.  Also,  as  temperature  and 
moisture  saturation  are  increased  (t*  -  0.5  and  250  °F)  the  value  for  the 
collapse  load  decreases.  The  overall  affect  of  the  material  degradation 
can  be  seen  more  clearly  in  Figs.  4.27  and  4.28.  These  figures  represent 
the  reduction  in  collapse  load  as  a  panel  with  no  cutout  is  altered  by 
both  a  cutout  and  then  temperature  arid  moisture  for  both  the  symmetric 
and  unsymmetric  cases  respectively.  As  Fig.  4.27  shows,  material 
degraded  panels  could  have  their  collapse  load  decreased  by  as  much  as 
&75  for  a  cutout  located  at  position  A  and  75.5%  for  a  cutout  located  at 
C.  In  other  words,  the  effect  of  placing  a  2"  x  2"  cutout  in  a  panel  and 
having  this  panel  subjected  to  high  temperature  and  symmetric  saturated 
moisture  conditions  is  to  develop  collapse  at  one-quarter  of  the  uncut, 
room  temperature  panel’s  buckling  load. 
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CONTOUR  LEVELS  ARE  IN  lOths  OF  MAXIMUM  DISPLACEMENTS 


.  4.24  Radial  Displacement  Contour  Plot  for  2 ‘  *  2*  Cutout 

at  Collapse  Load  of  198.5  Ibs/in  (0/ +45/-45/90)s ,  Cutout  A 
and  Symmetric  Moisture  Distribution,  t *  =  0 . 001 ,  and  Root 
Temperature  (80  *F) 
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Fig.  4.25  Radial  Displacement  Contour  Plot  for  2*  x  2"  Cutout 

at  Collapse  Load  of  169.0  lbs/in  (0/+45/-45/90)* ,  Cutout  A 
and  Symmetric  Moisture  Distribution,  t*=0.5,  and  High 
Temperature  (250  *F) 
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Fig.  4.26  Radial  Displacement  Contour  Plot  for  4’  x  4'  Cutout 

at  Collapse  Load  of  B2.25  lbs/in  (0/  +  45/-45/90) s ,  Cutout  A 
and  Symmetric  Moisture  Distribution,  t*=0.5,  and  High 
Temperature  (250  *F) 
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Table  4.5:  Material  Degradation  Results 


x  2"  Cutout 


Unsymmetric 

Moisture 


Symmetric 

Moisture 


Temperature 

( 0  F ) 

Nondimensional  Cutout 

Time  (t*)  Location 

Collapse 

Load 

Z 

Change 

Collapse 

Load 

Z 

Chanqe 

30 

0.001 

A 

199.4 

0.00 

193.5 

0.45 

C 

147.2 

0.00 

1 4  b  .  8 

0.2b 

NON 

514.4 

0.08 

513.3 

0.19 

30 

0.5 

A 

198.8 

0.30 

196.0 

1.71 

C 

147.0 

0.12 

145.3 

1.23 

NON 

510.7 

0.80 

506.8 

1.55 

250 

0.001 

A 

137.3 

6.09 

183.3 

7.32 

C 

138.9 

5 . 63 

1 3  b  .  3 

7.39 

NON 

4  79.2 

6.92 

4  7  6 . 6 

7.42 

250 

0.5 

A 

1 8b .  3 

6.57 

169.8 

14.84 

C 

137.9 

6.31 

126.3 

14.19 

NON 

452.6 

12.08 

432.6 

15.97 

t*  =  kt/h2 

Cutout  A 

=  199.4 

collapse  load  in  lbs/in 

Room 

temp . /no 

moisture 

Cutout  C 

=  147.2 

No  Cutout  :  514.3 
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Fi9.  4.28  Deduction  in  Collapse  Load  -  Unsr.oetric  noisture  Condition 
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V.  CONCLUSIONS 


The  following  conclusions  can  be  stated  based  on  the  results  developed 
from  this  study. 

1.  Collapse  characteristics  of  composite  panels  are  dependent  on  the 
ply  layup  and  the  size  of  the  cutout. 

2.  Two  ply  orientations  (0/  +  45/-45/90)s  and  ( 0/ -45  /  +  45/,:»0 ) s  were 
studied,  and  for  cutouts  located  along  the  vertical  axis  of  symmetry 
there  was  no  affect  on  the  collapse  load  due  to  the  two  ply 
orientations  (cutouts  A  and  8). 

3.  For  the  same  two  ply  orientations  stated  in  (2),  there  was  a 
reduction  in  the  collapse  load  of  17. bZ  for  a  2"  x2"  cutout  located 
near  the  side  support,  (eccentric  location). 

4.  Cutouts  moved  from  the  middle  to  the  side  can  produce  as  much  as  a 
2b. 7Z  reduction  in  the  collapse  load. 

5.  A  panel  with  a  small  cutout  (2"  x  2")  behaves  almost  linearly 
throughout  the  loading  and  at  collapse  produces  a  trough  like  that 


of  the  secondary  loading  path  (the  eigenvector). 


A  panel  with  the  larger  cutout  (4!l  x  4"),  due  to  the  larger  amount 
of  bending,  behaves  nonlinear ly  throughout  the  entire  loading  range. 

Material  degradation  effects  (high  temperature  and  saturated 
moisture  conditions)  can  result  iri  as  much  as  a  152  reduction  in 
collapse  load  for  a  panel  with  a  cutout. 

Material  degradation  effects  on  collapse  load  have  nearly  the  same 
percentage  change  no  matter  where  the  cutout  is  located. 

The  unsymmetric  moisture  condition  produces  up  to  a  t. .572  reduction 
in  the  collapse  load  of  a  panel  with  cutouts. 

The  symmetric  moisture  condition  produces  up  to  a  14.842  reduction 
in  the  collapse  load  of  a  panel  with  cutouts. 

A  panel  with  a  cutout  located  along  one  edge  subjected  to  saturated 
moisture  conditions  and  high  temperature  will  have  the  collapse  load 
reduced  by  as  much  as  75.52  when  compared  to  a  room  temperature,  no 


moisture,  uncut  panel. 


VI.  RECOMMENDATIONS  AND  SUGGESTIONS 


Panels  subjected  to  damage  (a  satellite  panel  striking  a  small 
object  in  space)  result  in  a  cutout  that  is  not  square  in 
nature  but  more  closely  resembles  a  circle.  Therefore,  work 
could  be  done  in  the  area  of  an  axially  loaded  composite  panel 
with  a  circular  cutout. 


(2)  Experimental  work  to  verify  the  results  of  this  thesis  would  be 
very  helpful  and  an  area  for  further  stu' 


More  work  should  be  done  with  the  STAGS  user -wr  i  t ten  subroutines. 
STAGS  has  the  option  for  the  user  to  enter  his  own  subroutines 
arid  this  would  allow  for  a  wider  range  of  problems  to  be 


investigated . 
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APPENDIX  A 


Appendix  A  consists  of  a  collection  of  plots  that  were  used 
in  the  analysis  of  the  information  presented  in  this  thesis 
These  plots  are  not  formally  discussed  but  are  included  as 
additional  information  for  the  reader. 
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Appendix  A1 


CUTOUT  LOCATION  A: 

Radial  displacement,  as  a  function  of  location,  taken  3,  6 
and  9  inches  from  the  loading  edge,  as  load  is  increased 


Moment  profiles  as  a  function  of  location,  taken  3,  6,  an 
9  inches  from  the  loading  edge,  as  load  is  increased. 


Fig.  Al.l  Radial  Displacement  as  a  Function  of  Location 
taken  3  inches  from  the  Loading  Edge 
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Fig.  A1.2  Radial  Displacement  as  a  Function  of  Location 
taken  6  inches  froe  the  Loading  Edge 


Triangle  :  61  1  b  s  / 1  n 

Diamond  ;  111  lbs/in 
Or  taaon  :  151  lbs/ 1  n 
f  Diamond  =  19?  lbs/in 

I 


Moment  Profiles  as  a  Function  of  Panel  Location 
taken  6  inches  from  the  Loading  Edge 
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Moment  about  X  Axis  (in-lbs/in  x  10'J) 


Triangle  :  61 
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Fig.  Ai.ll  Moment  Profiles  as  a  Function  of  Panel  Location 
faken  6  inches  from  the  Loading  Edge 
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LOADING  EDGE  DISPLACEMENT 
Cutout  location  A  -  A2. 
Cutout  location  C  -  A2. 


STUDY 
1  -  A2 
4  -  A2 


Cu tou t  loca t. ion  D 


A2.7  -  A2 


Fig.  A2.5  Loading  Edge  Displacement  vs  Applied  Load,  Cutout  C 
for  -45  Second  Ply 
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Fig.  A3. 3  W  Profiles  at  4.5  and  5.0  inches  fro*  the  Loading 
Edge  as  a  function  of  Location,  for  Cutouts  C 
and  +45  Second  Ply 


.00  12.00  24. 00  36.00  4<=  00 

Location  along  Panel  (°) 

Fig.  A3. 4  W  Profiles  at  4.5  and  5.0  inches  fro*  the  Loading 
Edge  as  a  function  of  Location,  for  Cutouts  C 
and  -45  Second  Ply 


Noaent  about  Y  Axis  (in-lbs/in  x  10'3) 


Fig.  A4.1  Hoaent  Profiles  at  x=4.750  and  x=7.250  inches  froa 
the  Loading  Edge  for  Cutout  A,  +45  Second  Ply,  at 
Collapse  Load  of  199.4  lbs/in. 
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Moment  about  X  Axis  (i 


Square  x=4.75  inches 
Star  x-7. 25  inches 


Location  along  Panel  (*) 


Fig.  A4.2  Hoaent  Profiles  at  *=4.750  and  *=7.250  inches  froa 
the  Loading  Edge  for  Cutout  A,  +45  Second  Ply,  at 
Collapse  Load  of  199.4  lbs/in. 
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Fig.  A4.3 


Hoeent  Profiles  at  x=4.750  and  x=7.250  inches  froe 
the  Loading  Edge  for  Cutout  A,  +4S  Second  Ply,  at 
Collapse  Load  of  199.4  lbs/in. 


Hoaent 


LOCATION  ALONG  PANEL 


Fig.  44. 6  Hoeen t  Profiles  at  x:4.750  and  x=7.250  inches  froe 
the  Loading  Edge  for  Cutout  C,  +45  Second  Ply,  at 
Collapse  Load  of  147.2  lbs/in. 


Location  along  Panel  (*) 


Fig.  A4.7  Moaent  Profiles  at  x=4.750  and  x=7.250  inches  froa 
the  Loading  Edge  for  Cutout  C,  -45  Second  Ply,  at 
Collapse  Load  of  178.6  lbs/in. 


Location  along  Panel  (•) 


Fig.  A4 . 8  Hoaent  Profiles  at  *=4.750  and  *=7.250  inches  froa 
the  Loading  Edge  for  Cutout  C,  -45  Second  Ply,  at 
Collapse  Load  of  178.6  lbs/in. 
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Fig.  A4 . 9  Moeent  Profiles  at  x=4.750  and  x=7.250  inches  frot 
the  Loading  Edge  for  Cutout  C,  -45  Second  Ply,  at 
Collapse  Load  of  178.6  Ibs/in. 


Fig.  44.10  Ho*ent  Profiles  at  x:4.750  and  x=7.25 
the  Loading  Edge  for  Cutout  C,  +45  Se 
a  Load  of  141  lbs/in. 


Fig.  A4 . 1 5  Moment  Profiles  at  x=4.750  and  x-7.250  inches  fro* 
the  Loading  Edge  for  Cutout  D,  +45  Second  Ply,  at 
Collapse  Load  of  178.5  lbs/in. 
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Fig.  A4.16  Noaent  Profiles  at  x=4.750  end  x=7.250  inches  froa 
the  Loading  Edge  for  Cutout  D,  -45  Second  Ply,  at 
Collapse  Load  of  147.2  lbs/in. 
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Fig.  A4.21  H*  at  the  Loading  Edge  of  the  Panel  at  Collapse 
Load  for  Cutout  C  and  +4S  Second  Ply 
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Fig.  A4.22  N*  at  the  Loading  Edge  of  the  Panel  at  Collapse 
Load  for  Cutout  C  and  -45  Second  Ply 
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